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Abstract 

We propose a model for local dynamics of a perturbed convex real-analytic Liouville- 
integrable Hamiltonian system near a resonance of multiplicity 1 + to, to > 0. Physically, the 
model represents a toroidal pendulum, coupled with a Liouville-integrable system of n non- 
linear rotators via a small analytic potential. The global bifurcation problem is set-up for the 
n-dimensional isotropic manifold, corresponding to a specific homoclinic orbit of the toroidal 
pendulum. The splitting of this manifold can be described by a scalar function on an n-torus, 
whose fcth Fourier coefficient satisfies the estimate 



O 



L-p\*M-\kw\ fcez n \{o}, 



where ui € M™ is a Diophantine rotation vector of the system of rotators; p € (0, ^) and a > are 
the analyticity parameters built into the model. The estimate, under suitable assumptions would 
generalize to a general multiple resonance normal form of a convex analytic Liouville integrable 
Hamiltonian system, perturbed by 0(e), in which case ojj ~ j = 1, . . . , n. 

1. Introduction and main result 

The main objective of this paper is to create a template to extend the theory for exponentially 
small separatrix splitting in Liouville near- integrable Hamiltonian systems near simple resonances, 
i.e. resonances of multiplicity 1, to the case of multiple resonances, of multiplicity 1 + m, m > 0. 
The interest in such a theory is dictated by the fact that the normal form theory and Nekhoro- 
shev theorem, resulting in exponentially long time stability (see e.g. is well developed for 

resonances of all multiplicities, whereas the exponentially small splitting phenomenon, resulting in 
similar exponents, has been quantitatively studied so far only in the special case of multiplicity one 
resonances. 

It is well known (see e.g. that a convex analytic Liouville near-integrable Hamiltonian 
system, with Hamiltonian 

H(p,q)=H (p)+eH 1 ( P ,q), (1.1) 

where (p,q) £ R n + 1 + m x T" +1+m are the action-angle variables on T*T n+1+m , where T = M/2vrZ, 
can be localized in the action space near a multiplicity 1 + to resonant action value p . Namely 
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suppose p is such that the kernel of the scalar product (DHq(p ), k), k G Z n+1+m is some 1 + m 
dimensional sublattice in x n+1+m . In this case, without loss of generality one can render DHq(p ) = 
(lo, 0) G M n+1+m . In addition, we further assume that lu G M. n is Diophantine, i.e. 

VfcGZ"'\{0}, \(k,Lu)\>$\k\- T ", (1.2) 

for some $ > and r„ > n — 1 (for n = 1 this obviously boils down to ^ 0). 

After a canonical change of variables, preserving the phase space bundle structure and time 
scaling the Hamiltonian (|1.1|) can be cast into the following normal form: 



H a t(p,q) = {^=,l) + ^(p,Q n{ p) + U(x ,...,x m ) + [/nf(q) + (p,g n{ (p,q))], (1.3) 

where p = (i, y , yi, ■ ■ ■ , y m ) G M n+1+m , q = ((p, x , x%, . . . x m ) G T n+1+m , Q n{ is a constant sym- 
metric matrix, and the pair (f n {,g n f) = 0(\/e) can be treated as a perturbation when e (suppressed 
in the latter notations) is small enough. Further in the paper, the bold typeface marks the n + l + m 
dimensional quantities. 

If one truncates the normal form Hamiltonian H n f by dropping the terms / n f and g ni in the 
formula (|1.3|) . the action l is flow-invariant. For l = 0, one can separate a natural system of 1 + m 
degrees of freedom, whose Hamiltonian can be written as 

K(y ,...y rn ) +U(x ,...x m ), (1.4) 

where K(y) is a symmetric positive definite quadratic form in y G R 1+m and U(x) - a scalar function 
on T 1+m . In the sequel, saying that some function is a "function on a torus" implies 27r-periodicity 
of this function in corresponding variables. 

In the simple resonance case m = 0, one can show that inherent in the normal form dynamics 
is the exponentially small separatrix splitting phenomenon, see |12| . |15j . 

In order to show how the exponentially small splitting theory can be built in the multiple 
resonance case m > 1, let us consider a simple model, which generalizes the so-called Thirring 
model for a simple resonance, see [5]. 

Namely, we study the following model Hamiltonian: 

1 n 

Hfj,(p, q) = (w, l) + - ^2 i| + Hi +m (y , ■ ■ ■ , y m , x , . . . , x m ) + nV(tp, x ,..., x m ), (1.5) 

where to is Diophantine, V is a real-analytic function on f n+1+m and /x is a small parameter. 
Specifically, for some strictly increasing sequence of positive reals Zo^i, • • • ^m, let H\ +m {y, x) have 
the following form: 

Hi +m (y,x) = K 1+m (y) + U 1+m (x), 

2 

K 1+m (yo, yi ,...,y m ) = lEi^of-' (1.6) 

Till 



Ui +m (xo,xi,...,x m ) = YT=o l i {U.T=i cosx j - 1 

Geometrically, the natural system (|1.6|) can be visualized as a "toroidal" pendulum, i.e. a particle 
of unit mass, confined to move on the surface of a "vertically standing" in M. 2+m torus of dimension 
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1 + TTi, with principal radii l m , . . . ,Iq, under the influence of gravity with the free fall acceleration 
equal to 1. Mechanically, the case m = 1 can be realized as a double pendulum, whose shorter arm 
of length Iq is attached to the terminal point of the longer arm of length l\ and moves in a circle, 
which rests upon the longer arm. 

On the energy level H^ m (0), the origin O = (0, 0) is a single hyperbolic fixed point, with the 
characteristic exponents 



1 

Suppose the sequence lj grows rapidly enough to ensure 



A > max(Ai, . . . , A m ). (1.8) 

In addition, we have to assume that 



inf 



3=1 



> 0, (1.9) 



where Z + denotes non-negative integers. 

The origin is connected to itself by a family of homoclinic orbits, in fact there exist homoclines 
representing each homotopy class on T 1+m for the geodesic flow generated by the corresponding 
Jacobi metric, degenerate at x = 0, see |3J. Some of these homoclinic orbits, or separatrices, 
are patently obvious: let xj = yj = 0, Vj 6 {0, . . . , m) \ {i}, Xi(t) = 4arctan e ±Xit . These orbits 
correspond to homoclinic geodesies forming the basis of the fundamental group of the torus T 1+m 
(modulo the sign in the exponential which bears witness to the reversibility of Hi +m , identifying 
the upper or lower separatrix branch, where yi retains its sign). Consider the orbit with i = 0, call 
it 7. This orbit leaves and arrives back to the fixed point in the maximum expansion/contraction 
direction, corresponding to the Lyapunov exponent Ao- In order to take both of the orbit's branches 
into account, let us represent 7 as follows: 

7 = {xi = . . . = x m = yi = . . . = y m = 0, y = 2 sin(x /2) = tp(x ), x <E (0, 2tt) U (2tt, An)}. 

(1.10) 

Observe that the existence of the two branches of 7, on each of which yo retains its sign, is reflected 
by 27r-antiperiodicity of the "separatrix function" ip: ip(xo) = — "0(^0 + 2ir). To reflect this fact, it 
will be further convenient to deal with xq G T2 = M/4ttZ rather than T = M/2ttZ. In particular, 
addition of values x is further meant to be mod {Air). 

Clearly, the orbit 7 belongs to both the unstable and the stable 1 + m dimensional invariant 
Lagrangian manifolds Wq S of the fixed point at the origin. If m > 1, the flow of the Hamiltonian 
H\ +m is non-integrable 1 . 

Global geometry of the manifolds Wq S is complicated. Locally near the origin however, the 
germs W^'w of the manifolds Wq S are diffeomorphic tom + 1 disks, tangent at O to the unstable 
and stable manifolds of the flow, linearized near the origin. 

1 The flow of -Hi+m should possess no global analytic first integral other than the energy, unless Ki+ m is diagonal 
and Ui+m separated, see |5j. Transversality of the intersection of the manifolds Wq' s along 7 (to be shown) is in 
turn an onset for non-integrability, see For the general variational approach to homoclinic trajectories in natural 
systems see [Sj- 
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We shall further show that 7 arises as a transverse intersection of the manifolds Wq' s . Let us 
call W"' s the localizations of these manifolds in the neighborhood of 7. As the orbit 7 takes off 
from/arrives at the fixed point in the maximum expansion/contraction direction, it itself turns out 
to be hyperbolic within the manifolds W^ ,s . Indeed, on the "vertical torus" in R 2+m , the coordinate 
directions the main curvature directions away from 7. 

Let us further change the notations (xo, yo) to (x, y), (xi, . . . , x m ) to z and (yi, . . . , y m ) to z and 
restrict \z\ < tq for some < To < 1. Then 

y 2 m 
Hi +m {y, z, x, z) = — 2 + /o(cos x-l) + ^2 

M o i=l 

The semicolon in the symbol 0^{z; xq) means that the term in question is 0(||z|| 4 ), uniformly in xq, 
|| • || standing further for the Euclidean norm, to be used intermittently with the sup-norm | • |. This 
notational convention will be used further on, the parameters following the semicolon often being 
omitted. 

Before formulating the main result, let us give some geometric description of what we are going 
to claim. Lifted into the phase space of the truncated Hamiltonian H^, where /x = 0, the orbit 
7 gives rise to an isotropic n + 1 dimensional invariant manifold, which is topologically a cylinder 
over the n-torus. Let us denote this cylinder as Co- Along Co, there intersects - degenerately in n 
directions corresponding to the rotators' variable <p — a pair of invariant Lagrangian manifolds W ' s , 
both containing an invariant whiskered n-torus To, located at (p, x,z) = (0,0,0). On the torus 
itself, the truncated flow is quasiperiodic, with the Diophantine frequency uj. Owing to the fact that 
Co has two branches, plus the fact that the trajectories on Co are bi-asymptotic to the invariant 
torus To, we shall technically refer to So as a "bi- infinite bi-cylinder" , yet tending to avoid this 
rhetoric, as much as possible. 

We study how the presence of the coupling term V in Q1.5|) is to affect the above described 
geometric structure and obtain qualitative estimates for the degeneracy removal effect. As far as 
the Hamiltonian H\ +m is concerned, the condition Q1.8JI results in local hyperbolicity of the orbit 
7 within the manifolds Wq S (recall that their localizations near 7 are denoted as W^' s ). I.e. the 
germs Wq^ oc will be contained in the closure of W^' s for the unstable/stable manifolds respectively. 
Let us denote W"' s = T n x W"' s , the lifting of the manifolds W^' s into the phase space of the 
truncated Hamiltonian H^, when ^ = 0. The manifolds TV"' S can be represented by their generating 
functions §>^' s (x, z) as graphs over the configurations space variables (<p, x, z), where (p £ T n , \z\ < r 
(for some small enough r to be determined) and x £ T2 \ (2tt — 5, 2n + 5) = [—2n + 5, 2ir — 5], for 
some positive 5 < 1. 

Then in the perturbed problem we are going to prove the existence of Lagrangian manifolds 
yju,s re p resen ting the analogs of the manifolds W"' s , as far as the Hamiltonian is concerned. 
Moreover, the homoclinic, or bi-infinite cylinder Co in the truncated system will give rise to a 
pair of semi-infinite cylinders Q u,s in the perturbed system, each perturbed cylinder containing 
the invariant whiskered torus T, the cylinders C" ,s themselves being contained in the Lagrangian 
manifolds W u,s respectively. The phase trajectories on W" will approach C" in negative time; in 
turn the trajectories on C M , in negative time will approach T at a faster rate. Similar orbit behavior 
will occur on W s and C s in positive time. This is the content of the structural stability theorem, 
Theorem 121 further in the paper. 

Moreover, the perturbed Lagrangian manifolds W u,s can be represented as graphs over the con- 
figuration space variables, by adding to the unperturbed generating functions S"' s (x, z) respectively, 



+ O i (z;x), (1.11) 
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some quantities S^' s (v?, x, z), which are both 0(h). Then let 

S u > s (<p, x, z) = §^(x, z) + S u /(<p, x, z) (1.12) 

denote the generating functions of the perturbed manifolds W u,s , respectively. 

In order to find these generating functions, we shall describe a series of canonical transformations, 
each of which explicitly takes advantage of the fact that the phase space is a cotangent bundle. In 
the sequel, any canonical transformation ^ will be determined by some automorphism a and closed 
one-form dS on the base space. I.e. all the canonical transformations dealt with herein have the 
following structure: 

¥ = *(a,5): { 9 u^'l-i , Ja , \ 

I V = (da) 1 p' + dS(q), v ' 

Observe that there is a natural semidirect product structure that on the pairs (a,S), induced by 
composition. 

Hyperbolicity of the orbit 7 does not suffice to prove Theorem |21 however: we also need a special 
non-resonance (yet not very restrictive) assumption (|1.9|) on the stability exponents of H\ +m at the 
origin, built into the choice of the arm lengths {lj}j=o,...,m- The latter assumption appears to be a 
very special case of the problem of analytic conjugacy between linearized and non-linear dynamics 
near a hyperbolic fixed point, see e.g. ^Hj> although for our purposes it suffices separating the 
dynamics in a single chosen direction only. 

The whiskered torus To and its local unstable and stable manifold Wq S 1oc are known to survive 
small perturbations without the assumption Q1.9|) . by the theorem of Graff, see jlUj . |18j . Namely, 
in the normal form (jl.3j) as long as u is Diophantine, U possesses a single non-degenerate absolute 
maximum, plus the upper left n x n minor of the matrix Q n f is nonzero, there exists a perturbed 
torus 7 where the flow is conjugate to the linear flow on its prototype in the truncated system. 

We further study the splitting of the unperturbed cylinder Co- In order to do so, we introduce 
the "splitting function" 

T>((p, x, z) = § u (<p, x, z) - § s (ip, x - 2vr, z), (1.14) 

which will be well defined for tp S T n , x £ [—2tt + 5, —5] U [5, 2ir — 5] (recall that addition of x is 
mod(4:Tr)), and \z\ < r. A critical point of D would yield a homoclinic connection to the torus 7, 
the gradient dT> being the "splitting distance" . 

As the manifolds W™' s for the Hamiltonian Hi +m intersect transversely at z = 0, the critical 
points of D will lie close to z = 0, and therefore, the magnitude of the splitting of the cylinder 
Co can be evaluated in terms of the derivatives D x ^T>(ip, x, z) at z = 0, in the properly adjusted 
coordinate chart 2 . 

The forthcoming Theorem ^ makes these claims precise. In order to formulate the theorem, let us 
introduce some notation and summarize what analyticity properties are required of the perturbation 

V((p,x,z). 

For real r, a > and j = 1, 2, . . . (j = 1 usually being omitted) let 

Mi ^ {C € : HCII < r}, 

it d ^ f {( e : 5ftC G T-?, |3fC| < a}. 

2 This is the only instant in the argument of this paper, where the built into the model transversality of the 
intersection of the manifolds W^' a comes into play. 
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For x G T2, define a conformal map s and some associated quantities as follows: 

f x d( 

s{x) = J x(s) = ip[x(s)], e = yx{s). (1.15) 

Recall, in the model studied ip(x) = 2sin(x/2). The map s(x) takes (0, 4tt) to R U M + «7r, and 
the change (x,y) — ► (s,e) is canonical. The function x(s) is 2-7ri-periodic and has singularities at 
S = ±fi 

Fix some To 3> 1. By construction of the map s, for any and p £ (0, 7r/2) any T 6 [To/2, To], 
the quantities a;(s),x(s) are holomorphic functions in the set tlr,p C C/2iri, obtained by throwing 
out of C horizontal rectangles with half-axes (2To — T) x (vr/2 — p), centered at ±?i. Namely, let 

nT,p = n TjP u-n TiP , 

n T , P = {3fo<T, |3fs| <p}U{Ks<T, |3fs-7r| < /?} U {5fts < T - 2T }; also let (1.16) 

n Tp = { s G C : \Us\ < T,\Ss\ < p}. 

The domains Ht,p are further referred to as semi- infinite bi-strips, their size increasing with (T, p), 
with p < \. Bi-strips t^T,p are bi-infinite, while IIt.p is simply an origin-centered horizontal rectangle 
in C, with semi- axes (T,p). 
Let 

e CT ,T,p = T™ X n T , p , C,T,p,r = C^T.p X IB™ (1.17) 

(and in the same fashion C a ,T,p, &a,T,p,r or Q a ,T,p, <£-a,T,p,r) be referred to as complex semi-infinite 
(bi-infinite or finite) bi-cylinders for C and extended bi-cylinders for the notations £. In qualitative 
argument, the analyticity indices as well as "bi-" rhetoric are avoided. 
Let us now quote the main assumption. 

Assumption 1 Assume the non-resonance conditions M.°J\) and Kl.ty . Suppose the real-analytic 
function V((p,x,z) is such that V(<p,x(s),z) is holomorphic and uniformly bounded by 1 in 
£a ,T ,po,ra f or some initial set (o"o, To, po, ro) of analyticity parameters. 

The main result of the paper is the following theorem. 

Theorem 1 Under Assumption^ take T = To — 1 and any positive p < po, a < ctq, let 6 ~ logT. 
Suppose 

r < Cl mm[(p - p),(a -a)\, p < c 2 [r#\u\- 1 (a - a) T "] 2 , (1.18) 

for some constants c\% > 0, determined by the separatrix function ip as well as the quantities 
n, r„, m, cj , T , po,r , lo,... l m . 

i. Some level set of H^, with energy O(p), contains an invariant partially hyperbolic n-torus 
7, where the flow is conjugate to linear, with the rotation vector to. At the torus 7, there 
intersects a pair of isotropic manifolds C u and C s , which are contained respectively in the global 
unstable and stable manifolds of 7. The manifolds Q u and Q s are contained respectively in a 
pair of Lagrangian manifolds W u ' s , which are graphs of closed one-forms, with the generating 
functions § u ' s (ip, x, z), as in 11.12\) respectively, such that the quantities Sp s ((p, x(s), z) are 
holomorphic and uniformly bounded by 0{p) for ((p,x,z) S <£ a ,T,p,r- The cohomology classes 
£u,s g H l (T n , R) = W 1 of the one-forms dS u ' s are equal to each other. 
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ii. The distance between the manifolds W u ' s can be measured by the exact one-form dD, defined 
by {1.1$ . There exist a coordinate chart ((p',x',z') G T n x [8, 2ix — 5] X BJ?, obtained by a 
near-identity change of variables (ip',x',z') = a(<p,x,z) from the original coordinates in hl.5\) . 
such that in the chart ((p\x\z'), the function T> satisfies the following PDE: 

*W^ + (»,j£) + V,L'V>])=0, (1.19) 

where U is a linear first order differentiation operator. 

Hi. In the above chart, namely for (ip',x',z > ) £ <£-o,T,p,r, the function 2) (</?', x', z') is bounded by 
0(/j,). Let 2) (<£>', x 1 , z') = T>o(x',(p') + 0(z'). Then the quantity T)q (x',ip') can be written as a 
2tt -periodic function on T n of 

a = (f/ - us(x'), i.e. T> (x',ip') = &(a). (1.20) 
iv. The manifolds W u and W s intersect at least 2n + 2 orbits, biasymptotic to 7. 

Let us show for the moment that the conclusions (iii) and (iv) of Theorem ^ are straightforward 
consequences of (i) and (ii). Indeed, (|1.19j) implies that Do(x',(p') has to satisfy the linear PDE 

d2) , dV 

^ (x) a7 + ( "'a7 ) = ' 

and ip(x')jjp- = Jjj where s(x') comes from (|1.15j> . Then T>o(x',(p') = 6(a) follows, as the form 
dT>Q is exact, i.e. T)q(x' , ip') is 27r-periodic in ip' £ T n . It follows that the set of the critical points 
of the function &(a) in the coordinate plane z 1 = determines the trajectories, biasymptotic to 
the torus T. The minimum number n + 1 of critical points of & is the Ljusternik-Schnirelmann 
characteristic of T n , which equals n + 1. It gets doubled in the statement (iv), because one can 
restrict x £ [—2tt + 5, —5], considering the lower separatrix branch and replicate the statement (ii). 

Theorem ^ has an immediate corollary, implying the estimate claimed in the Abstract and 
exponential smallness of the splitting distance if uj — * for a small e. We do not elaborate 
on various parameter relations in this paper (they can be all made to depend on e in order to 
approach the lower bounds' problem for the exponentially small splitting distance) as the situation 
here would be the same as it is in the simple resonance case, regarding which see e.g. [1] and 
the references therein. The willing reader may synthesize these relations using the fact that the 
parameter relations in the forthcoming technical statement Theorem are supposedly optimal. 
Going carefully through the latter theorem, one can derive what precisely the symbols 0(fj,) imply 
regarding the other parameters of Theorem ^ 

Corollary 1.1 For k £ Z™ \ {0} ; the Fourier coefficients &k for of the function 6(a) satisfy the 
estimate 

|6fc| < O(n) ■ e -p[M-[*[» (1.21) 

If (<p,x,z) are the original coordinates and ((p',x',z') = a'(ip,x,z) is the change of variables, de- 
scribed by Theorem^ (ii), then for all (<p,x,z) £ T n x [5, 2ir — 5] x [—r,r] m , one has a uniform 
bound 

\V oa'(ip,x,z)\ < 0(h)- e-"^-^. 

fceZ"\{o} 
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Observe that Theorem ^ is also valid in the simple resonance case m = 0, when the manifolds S and 
Q u ' s are Lagrangian, rather than isotropic. The simple resonance case has been exposed in detail 
in |15j . see also the references contained therein. In order to keep the ideas clear and not to rival 
the latter reference length-wise, intermediate steps in deriving many estimates in this paper are 
omitted, no intermediate analyticity parameters are explicitly introduced, and the constants c\ 2 
can get smaller from one statement to another. The reader is often referred to ^H^, as many aspects 
of the routine here mimic the m = case. The proof of how Corollary 11.11 follows form Theorem ^ 
can also be found in the latter reference, as well as JJ], |12| . 

Remark. In the perturbation V((p, xq, . . . , x m ) in (|1.5[) . let Vo(<p, xq) = V((p, xq, . . . , 0), and suppose 
it vanishes if xq = 0. Then one can formally set up the Melnikov integral 

/+00 
V (a + ujt,x (t))dt, aeT", (1.22) 
-00 

and ask whether this quantity is a bona fide C 2 approximation for 6(a). (A routine calculation 
shows that setting to — ► results in the Fourier coefficients Wl^, equal to the right-hand side of 
the bound p.21j) . with p = ~.) We do not study this undoubtedly important issue here (see. e.g. 
[I], ^2] f° r thorough discussion). However, the theory developed further suggests that there are 
no extra difficulties arising in this respect in the multiple versus the simple resonance case. In 
particular, the "easy" case n = 1 involving no small divisors should be similar in this respect to the 
Thirring model, cf. [S]. 

Remark. Note that contrary to the simple resonance case, where there exists a large body of 
literature on exponentially small bounds for the splitting of separatrices, see e.g. ^2] and the 
references therein, multiple resonances have been usually approached via the normal forms, alias 
averaging method. The latter technique (see ^2] , |14j ) is not very explicit geometrically, however 
as ^1] points out, it does enable one to obtain exponentially small upper estimates with sharp 
constants, which come from dynamical considerations regarding the analyticity domains, if not to 
relate these estimates directly to the splitting of separatrices. 



2. Unperturbed system analysis 

In this section, for the sake of clarity, we confine ourselves to the case m = 1 only; the extension to 
m > 1 is transparent. Thus, in this section, let Iq = 1, l\ = I > 1. Let us further assume that I is 
such that A = ^±1 < 1 anc i f or k £ Z + 

\k\-l\ >—, Vk g Z+, (2.1) 

a particular case of the condition (|1.9jl . 

The Hamiltonian Hi +m given by ()1.11|) . for m = 1 turns into 

y 2 z 2 z 2 

H 2 (y,z,x,z) = — + (cos 2 - 1) + ^2 ~ ( l + cosx)y + 4 {z\x). (2.2) 

Let us further change y — > ±ip(x) + y, recall that ip{x) = 2sin(x/2). Clearly, the choice of the 
sign as + corresponds to localization near the orbit 7 as part of the unstable manifold of the fixed 
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point O, while the — sign would imply doing it near 7 as part of the stable manifold. Let us call 
the resulting Hamiltonians i?2,±V> as follows: 

xp 1 z 2 

H 2 ,±i>(y,z,x,z) = ±yip(x) + Y + 2P ~ ^ + cosx )y + Oi(z;x). (2.3) 

Observe that the Hamiltonians H 2t ±ih nave resulted from H 2 after the canonical changes with the 
generating functions 

S± = ±J%(C)d(, (2.4) 

with the + sign for the stable and the — sign for the unstable manifolds, respectively. Further 
calculations will be quoted for mostly H^ +^p = H 2 ^ only. 

The Hamiltonian H 2 ^ is now a function of x £ T2 = R/47rZ, rather than T, with two singular 
points, where x = 0, 2tt. To identify them H 2 ip, retains a symmetry: 

H 2 tf(y, z, x, z) = H 2 ^(y + 2ip(x), z,x + 2ir, z). (2.5) 

This symmetry was called the "sputnik property" in 15 and was used to validate the analogue of 
the claim £ u = £ s of Theorem ^ in the case of a simple resonance. Here, as we are dealing with 
the specific model, we will not use this property explicitly to prove this claim, but rather will later 
compare the pair H 2 ^ (which is anyway tantamount to the same trick, used in the beginning of 
section 4). 

To identify the unstable/stable manifolds W^' s of the orbit 7, we will be looking along the 
x-axis at the unstable/stable manifold of the singular point at x = for the Hamiltonians H 2t ±^p, 
respectively. This proves convenient and suggests that in general the regular description of the 
manifolds Wj' s is likely to fail in the neighborhood of x = 2tt. 

Let us linearize the flow of the Hamiltonian H 2 ^ near the orbit 7, whereupon x(t) = xo(t) = 
±4arctane*. For the infinitesimal increments (x,y,z,z), one gets the system of equations 

x = D^[x (t)]x + y, y = -D^[x (t)]y, (2.6) 



k = z/l 2 , I = (l - 1 + tanh 2 t)z. (2.7) 

The tangent space to W" at the points on 7 will be spanned by the vectors - solutions of the latter 
system of linear ODEs, which vanish as t — > —00. 

The two pairs of equations ()2.6|) . (|2.7f) are uncoupled. As far as (|2.6|) is concerned, there is 
an obvious solution x(t) = xq{€) ~ 1/ cosht, y(t) = 0, which vanishes at both t — * =Foo. I.e. one 
tangent direction to at a point (x,z,y,z) = (x, 0,0,0) is always (1,0,0,0), in the direction 
collinear with 7 itself. 

Equations (|2.7j) will clearly have no solutions vanishing at both i±oo, as the coefficients therein 
retain their sign for all t. However the system certainly does have a solution (z u (t), z u (t)), defined 
for t < Tq for some To S> 1, which as t — > —00 approaches the trivial (z,z) = (0,0) (as well as 
another solution, which is defined for t > —To and vanishes at t — * +00). 

To construct the unstable solution (z u (t) , z u (t)) , one may set ^ = ip{x)^, 1 — tanh 2 t = cosx 
and construct the germ of the solution in question locally as a Taylor series in x near x = 0; by 
linearity of the equations <|2.7|) and boundedness of their coefficients, the continuation of these germs 
over a finite time interval does not itself pose any problem 3 . 

3 As a matter of fact, equations l|2.7^ represent a second order linear ODE of generalized Legendre type, whose 
general solution can be found explicitly in terms of the associated Legendre functions. 
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Observe that given (z u (t) , z u (t)) , one can let {z s (t), z s (t)) = (z u (—t) — 2ir, —z u {—t)) for the result 
of the similar procedure with respect to the Hamiltonian H 2 — ,/>■ For no t E [—To, To] can the vectors 
(z u (t), z u (t)) and (z s (t) + 2tt, z s (t)) be parallel, or there would exist a solution of ()2.7|) . biasymptotic 
to zero. For the Hamiltonian H 2 from (|2.2|) . the existence and transversality of the intersection 
along the orbit 7 of a pair of manifolds W^' s (defined in the neighborhood of 7) essentially follow. 
A quantitative statement of this fact is to be given shortly. So far observe by comparing the 
coefficients in the linear equations (|2.7|) retain the sign, both vectors (z u (t) , z u (t)) and (z s (t), z s (t)) 
in the (z, z) plane never have a slope too close to horizontal. More precisely, their slope in absolute 
value will be contained in the interval [ly/l — 1, lyl + 1], which can be seen from (|2.1U|) below. 

Let us show how the manifold can be constructed, the analysis for gets modified in the 
obvious way. Make a change z — ► z + X u (x)z, where X u (x) determines the direction of the solution 
vector, vanishing at x = (i.e. t — > —00). The quantity X u (x) G [ly/l — 1, ly/l + 1] is well defined 
for x G T2 \ (2tt — 5, 2tt + 5) for some < 5 < 1, where <5 « In To. 

The change z — > z + A u (x)z is not canonical, to make up for it one also has to change y — ► 
y + ^—^^-z 2 - In other words, this is a canonical change with the generating function 

§l (x,z) = ^X u (x)z 2 . (2.8) 
Then the Hamiltonian i?2,^ in <|2.3|) transforms to 

H 2 , u (x, y) = y^{x) H h -A u (x)z 2 + r 2 X u (x)zz + y0 2 (^; ar) + 4 (z; x), (2.9) 

where 

dX (t\ 

X(x) = i>{x)— |p - (Z + cosx) +/- 2 A2(a;). (2.10) 

It follows by construction - or directly from (|2,7|) - that X u {x) = 0. Thus the quantity A u (x) = 
l~ 2 X u (x) multiplying zz in (|2.9() is always positive, never exceeding A = A u (0); recall that A < 1, 
also cf. (EUD - 

Finally, the last two terms in ()2.9|) can be regarded as a perturbation, provided that r is small 
enough. 

The phase space of the Hamiltonian H2 )U is T*((T 2 \ (2tt — 5, 2tt + 5)) x [— r, r]). If there were 
no two last terms in (|2.9|) . the manifold W" would be given by the zero section (y,z) = (0,0) of 
the bundle. However, for small r, the last two terms in (|2.9|) can be regarded as a perturbation and 
dispensed with, owing to the following lemma. 

Lemma 2.1 Given p < p$, T < Tq — 1, there exists a constant c\ > 0, depending only on the 
•parameter set (po, To, ro) and X, such that for r < c\(po — p), there exists some reals 5,n = 0(1) in 
(0, 1) and a canonical near-identity transformation VP", such that the Hamiltonian \2. ,9J) can be cast 
into the following normal form: 

H 7>u (y, z, x, z) = yip(x) + A u (x)zz + 2 (y, z), (2.11) 

valid for \y\, \z\ < k, \z\ < r and x such that $tx £ [—2ir + 5, 2ir — 5] and s(x) G ^T,p- 

The transformation for p = (y, z) and q = (x, z) can be written in the following form: 

n = * u M,K)-- [l = i [id + Ju {qf)] - y + (2-12) 
where the quantities b^(x,z) and §^(x,z) are both 02(\x\ + \z\). 
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We omit the proof of the Lemma, as it follows as a particular case of the forthcoming Theorem [2J 
in the case when there are no ^-dependencies. The smallness condition r < c\(po — p) follows after 
routine, but careful examination of the proof of Theorem see also the quantitative estimates in 
Theorem [21 As a matter of fact, Lemma 12. II still holds if the term 0±(z; x) in (|2,9j) gets replaced 
by 3 (z;x). 

Observe that repeating the argument for the Hamiltonian the latter would be cast in the 

following form: 

H^ :S (y, z, x, z) = -yip(x) - A s (x)zz + 2 (y, z), (2.13) 

where A s (x) > and equals A at x = 0. In order to get (|2.13j) . the analog of Lemma 12.11 would 
be preceded by a canonical transformation with the generating function §i f0 (x,z) = ^X s (x)z 2 , cf. 

Let us summarize the results of the analysis in this section by the following proposition 

Proposition 2.2 Forr small enough, the unstable/ stable manifolds W^' s for the Hamiltonian h2.2]) 
can be represented as graphs over the variables (x,z), via the generating functions 

S^ = S+'-+S^ + S^, (2.14) 

respectively, the representation being valid for $tx G T2 \ (2ir — 5, 2tt + 5), s(x) £ ^p,T, S ~ log T and 
z G B™. Both §"' s vanish to the second order at (x,z) = (0,0). The intersection of the manifolds 
W^' s along the orbit 7 is transverse, for x £ [5, 2tt — 5] U [2tt + 5, 2tt — 5] . 

The Hamiltonian can be cast into the forms 12. lip . \2.13\) via canonical changes ^"' s 

respectively, where *$>™' s = ^' s (ar' s , S"' s ), with the near-identity diffeomorphisms a"' s = id + 6" ' s 
such that br' s both vanish to the second order at (x,z) = (0,0). 

The next and most important step towards proving Theorem^is to get the generating functions 
S^' s in ()1.12j) by developing the structural stability theory for a class of Hamiltonians, which would 
include 

1 n 

H u>s = («,i> + -2 l ) + H 7,- + Vu,s(f, x, z), (2.15) 
1=1 

with • = u,s, respectively, and the perturbations V U)S satisfy Assumption ^ We can also assume 
that we can estimate the partial derivatives of D XyZ V UyS at (x,z) = (0,0) in terms of the norm of 
V u ,s, as V U}S itself should be analytic for |x| < 5 and \z\ < ro- 

3. Structural stability theory 

The classical structural stability theorem in problems with small divisors is due to Kolmogorov 
|llj . alias KAM. It establishes stability of geometric objects - invariant Lagrangian tori with quasi- 
periodic flow thereupon - with respect to perturbations of non-degenerate Hamiltonians localized in 
the neighborhoods of these geometric objects. If one is after other geometric objects, say whiskered 
tori of lower dimension, with quasi-periodic flow, there is a theorem of Graff, JHJ- It was proposed in 
|15j that a proper geometric object to look at in order to set up the splitting problem near a simple 
resonance is a semi-infinite cylinder over a torus. In this section the structural stability theorem 
from is given extra development. Namely, we deal with "extended cylinders", see ()1.17j) . which 
appear to be the proper geometric objects to study in order to describe the separatrix splitting 
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at multiple resonances. Naturally, all the above mentioned theorems, as well as the forthcoming 
Theorem|21 which is the main result in this section, are in the realm of the abstract implicit function 
theorem framework of Zehnder ^S], |2L)| . 

As this section is the most technical one, the notations in it (hence also in section 5 and the Ap- 
pendix) are largely self-contained. We study the following Hamiltonian, with (p, q) = (t, y, z, p, x, z) 
(to justify the earlier made claim that it suffices to consider m = 1, let us now take z = (zi, . . . , z m ), 
for any m): 

H u (p, q) = \o^{x)y + (u, i) + (z, A(x)z) + 2 (p; q), (3.1) 
under the following basic assumptions: 

i. uj is Diophantine, for all p, the matrix D^ L H W (Q, 0, 0, p, 0, 0) is non-degenerate. 

ii. For all x E T2 \ (2tt — 5, 2ir + 5), the real parts of the eigenvalues of the diagonalizable matrix 
A(x) lie in the interval (0, Ao); moreover A(0) = diag{\\, . . . , A m ) and the condition (jl,9[) is 
satisfied. 

Clearly, the simultaneous change (Ao, A) to (— Ao, —A) is not going to violate the principal conclu- 
sions of this section. 

Notation 

Technically it proves convenient to deal with the non-compact "energy-time" coordinates (e,s), 
introduced by ()1.15j) . rather than the coordinates (y,x); some notation and formalism are being set 
up further. 

Let QSo-ClP) be the Banach space of bounded 2-7r-periodic scalar functions in each variable, real 
analytic in T£, with the sup- norm. 

Let x E T2. Rather than dealing with a fixed ip{x) = 2 sin(x/2), let us introduce it axiomatically, 
as a real-analytic function, such that V(0) = 0, Dip(0) = 1. Suppose ip(x + 2ir) = —ip(x) and ~4>(x) 
has no other zeroes on the real line, but integer multiples of 2ir. 

Define a conformal map s(x) via ()1.15j) . The map s(x) takes (0, 4ir) to MUM + i7r and the change 
(x, y) — > (s, e) is canonical. 

By construction of the map s, there exists some 3> 1 and p E (0, vr/2) such that for any 
T G [T^/2,T^j] the quantities x(s),x(s) are holomorphic functions in the set IIt,p C C/2ni, defined 
by ()1.16j) . In addition, and this is possible by the properties of ip(x), let us suppose that p is such 
that for any s E II^p, there exists a pair of constants c^, such that 

c^e-l s l < \ X (s)\ < C^ e -I S L (3.2) 

To deal with Hamiltonian functions in T*(£ let us introduce some function spaces. For more details, 
see [T5j . 

The function spaces will be characterized in terms of the analyticity parameters, accommo- 
dated into parameter vectors p as follows. Let p = (a,T,p) E Introduce partial order 
p' = (a', T', p') < p if a' < a,T' < T, pf < p. If p' < p and |p - p'| = inf (a -a',T- T',p - p') > 0, 
write p' < p. Addition of parameter vectors, as well as multiplication by positive reals is de- 
fined component- wise, as well as the difference p — p' for p' < p. For A E (0, |p|), the notation 
p' = p — A means that A has been subtracted from each component of p. In the sequel the com- 
ponents and dimension of the parameter vectors p may vary; more often than ever we will have 
p = (<r,T,p,r)eR% + . 
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Let 58£(C) be the Banach space - with the sup-norm - of bounded holomorphic functions u 
on Cct,t,p, such that u(ip,s) = u[ip,s(x)] = u((p,x), where u is bounded and holomorphic in the 
Cartesian product of T CT and the pre-image of the strip H-t,p under the map s(x), and u(p, x) 
vanishes to the jth order at x = (the index j = being omitted). Component- wise sup- norm | • | p 
or the equivalent Euclidean norm || • || p is used for vector functions. 

If u((p,s) G *Bp(C), a multiplier x J ( s ) can De factored out, i.e. 

ufas) = x ] {s)v(ip,s), vG<B p (C), \v\ p « \u\ p , (3.3) 

with constants depending on the fixed quantity x being henceforth absorbed into the symbols <, ~, 
see also (|3.2j) . For u G 23p(C), there exists a unique decomposition 

ufas) = u fa+ Ul fas), where u G *B a (T n ), u\ G 33p(C). (3.4) 

Using it, define the average {u) "at infinity" as 

(u) d = / u (tp)dip. (3.5) 

For u G 5Sp(C), there is an estimate: 

l«(y>*)| ^ eS klp- 

Let us also introduce a function space ?Bp(C) = *B (T (T Tl ) 
infinity as follows: 

ufa s) e 55^(6) iff ufa s) = vfa s) G <B P (6). (3.7) 

The norm on 25p(C) is defined as \v\ p . By (|3.4|) and (|3.7j) . for -u G 5Sp(C) there is a decomposition 

u(y>,s) = i>ofa)/x(s) + vi fas), for some vi G » p (e). Also let 53j n ' A) (e) = [<B P (e)f x 58^(6). An 
element of this space describes a vector field on C as well as a map a((p,s) of S p into S p /, (with 

p < p'for the map to be well defined). Namely, if g G Q5p n ' A \c) is a vector field and a is such a map, 
then the "new" vector field da~ l g o a is in 93^ (C), see ^Sj. It is legitimate to use the Cauchy 
formula to estimate partial derivatives of u G Q3 A (6), i.e. \du\pi < A _1 |u| p , where A = p — p'. 

All the above notations extend in an obvious way to functions on £, by adding a component r 
to the parameter vector p, and considering absolutely convergent Taylor series in z G Wf 1 with the 
coefficients in the corresponding spaces of functions on C. For instance 53p(<£) (with p = (a, T, p, r)) 
becomes an extension of the space 25p(C) (with p = (a,T,p)), and the notation Q$v rl ' A ' m ) (£) extends 

25p A ^(C). The norm in the extended spaces, such as 93p((£) is the sum of the Taylor series in z, 
where the moduli suprema have been taken for all the coefficients. The quantity r will not appear 
explicitly in the estimates in this section, getting absorbed in the < symbols, e.g. for u G 25 p ((£), 
we have |L> 2 u| 2= o < \u\p- To bring it to terms with the fact that r in Theorem ^ is actually quite 
small, cf. (|1.18|) . observe that r would further come into play only when the functions' derivatives 
are evaluated at z = via the Cauchy inequality. But for the functions Theorem ^ is dealing with, 
these derivatives can be estimated in terms of the quantity ro, which is 0(1). The same should 
be said about the parameters (<r, T, p) which are all supposed to be independent of the parameter 
characterizing the perturbation size. 



(3.6) 

x 5B P (C) of functions unbounded at 
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The notation (u) for u G *Bp(£) implies that z has been set to zero, cf. (|3.5|) . Thus, for 
u G 95j(£), there is a uniform estimate, cf. (|3.6j) : 

K^s,*)! i$ e s |u| p . (3.8) 

Hamiltonian functions on T* (£ are given by absolutely convergent Taylor series with coefficients in 
Q3p(£),inp= (i,y,z) = (i, e/x(s), 2), inside a complex ball B™ +1+m . Notation-wise (e, s) = (0, — oo) 
corresponds to y = 0. Let 2J K) p(T*(£) be the space of such Hamiltonians, the norm being the sum 
of the Taylor series in p, where the norms have been taken for all the coefficients. 

Structural stability theorem 

What follows is a non-technical formulation of the theorem to keep its content transparent. 
Theorem 2 Consider the following Hamiltonian H w G 23 Kj p(T*£): 

H u {l, e, z; cp, s, z) = const. + A e + (z, A(s)z) + (u, l) + 2 (p; q), (3.9) 

with p = (i,y/x(s),z). Assume the following: 

i. uj G M. n is Diophantine and the matrix D^H u (0, 0, 0; cp, — oo, 0), is non- degenerate. 

ii. Ao > 0, i/ie rea/ parts o/ a// the eigenvalues of A(s) G ?8«(C), Vs Zie m £/ie interval (0, Ao); 

ra. in i/ie decomposition A(s) = Ao + Ai(s), with Ai(s) G QSj(C), one has Ao = diag(\i, . . . ,X m ), 
with < 5RAj < Ao, Vj = 1, . . . , m, and the condition M.9\) is satisfied by {Ao, . . . , A m }. 

Then H u is structurally stable, via a canonical transformation 

and for any p' < p, the transformation a = id + 6, with b G ^8jj' A (C). The one-form dS is defined 
by the generating function S(q) = (£, tp) + So(<p, s), with £ G R n , So G *8 p /(C). 

Let 

H = H " + V > (311) 

be a small perturbation of the Hamiltonian (|3.9j) . In the perturbation, suppose / G 53p((£) and 
fif G 5Sp n ' A,m ^(£). How small the perturbation should be is stated in the forthcoming technical 
version of Theorem [21 Theorem |21 

Remark. An important consequence of the analytic set-up to be used further is local uniqueness. 
I.e. given the pair (H^,V), the pair (a,S) defining the conjugacy transformation ^ in ()3.1(Jj) is 
unique. 

The proof of Theorem |2] is given in section 5. Let us now discuss some implications of the theorem, 
in the coordinates (cp, x, z), cf. (|1.15jl . where cp G T n , x G T2 \ (2ir — 5, 2-it + S), \z\ < r. Let S, So be 
the expressions for the generating functions S, So from Theorem in these coordinates. 
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Corollary 2.1 The Hamiltonian H in h'J.ll\) . as a function of(i,y,z,(p,x,z), on some energy level, 
possesses an invariant Lagrangian manifold <£, given by the graph of the closed one-form dS, where 

S(ip,x,z) = (£,¥>) + S Q (<p,x,z), (3.12) 

and So is 2ir -periodic in the variable p. The manifold <£ contains a partially hyperbolic invariant 
torus 7, which in turn is contained in an invariant cylinder C = T n x [—2ir + 5,2ir — 5]. 

If the perturbation (f((p,x(s),z)),g(ip,x(s),z)) in \3.11\) is such that f = C>2(\x\ + \z\) and 
g = 0\{\x\ + \z\), then £ = and the energy value on the manifold £ coincides with the value of the 
unperturbed Hamiltonian thereon. 

Indeed, the first claim follows from (j3.9j) and (|3.10j) by setting in the latter formula p' = 0. Fur- 
thermore, if S(ip, x(s), z) = S(ip,s,z), where the latter comes from Theorem^ then the invariant 
cylinder C arises by letting z' = in S'(tp' ', s' , z') = S o a^ 1 ^' , s' , z'), where the transformation a 
also comes from Theorem [21 The torus 7 arises by further setting s' = —oo. 

The second claim follows by observing that a special perturbation, as described in the Corollary, 
does not affect the invariant torus at (p, x,z) =0 and local uniqueness. Alternatively, one can 
verify this claim by carefully inspecting the proof of Theorem 

4. Conclusion of the proof of Theorem [l] 

Combining the claims of Proposition 12.21 Theorem [21 as well as Corollary 12 . 1 1 applied to the Hamil- 
tonians H U)S in (|2.15j) . one immediately establishes the claim (i) of Theorem ^ but for the fact that 
the invariant manifolds W u ' s lie on the same energy level and the fact of equality of the cohomol- 
ogy class representatives G W 1 = H l (T n ,M). Both facts however easily follow by observing 
that all the generating functions in (|2.14|) vanish to the second order at (x, z) = (0,0), where the 
unperturbed invariant torus is located, so one can use Corollary 12.11 

Namely, let a"' s extend the diffeomorphisms a"' s in Proposition l2.21 acting as the identity on the 
(^-variables; let ^' s (a^' s , §>™' s ), be the corresponding canonical transformations. Let ^^' s (a^' s , S^' s ) 
be supplied by Theorem [21 being applied to the Hamiltonians (|2.15[) . where the quantities a^' s ,S^' s 
are viewed as the functions of (ip, x, z) rather than (ip, s, z). Let H%' s be the results of conjugacy: 

TTU,S TT tyU,S q,U,S 

respectively for the unstable and the stable manifolds. 
Consider the Hamiltonian 

H' = H%o (l^ 1 o 

By the properties of the pair (gv' s , S"' s ) described by Proposition 12.21 it follows that H' = H'^ + V, 
where V = (/' ,g') is such that / vanishes to the second order and g' to the first order at (x, z) = 
(0,0), while H'u can be regarded as unperturbed Hamiltonian, in the sense of Theorem [21 This 
implies that by Theorem[2]and Corollary 12. 1[ there exists a transformation ^f'(a', S'), which nullifies 
the perturbation V', and the one-form dS' is exact, i.e the corresponding £' = 0. 
Thus 

H'oV'= [H^ o *«) o [(tt-)- 1 o o tf« o (tt*)" 1 o o ^ 

i.e., by uniqueness, the application of Theorem[2]to the "stable manifold" Hamiltonian H s = H^oty^ 
is effected via the canonical transformation 
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This transformation is still of the form (|1.13|) . Besides, the corresponding generating function will 
contain a single "non-exact" term (£, u ,ip), supplied by VP" as (in the sense of the template 1(1.13(1 1 
the rest of the transformations in the above chain are effected by exact one- forms. This proves the 
claim (i) of Theorem ^ 

To prove the claim (ii) of the theorem, substitute p = dS u (q) and p = dS s (q) into the Hamil- 
tonian ((1.5(1 . subtract the result of the latter substitution from the result of the former one. After 
substraction has been done, all the momentum-independent terms are gone, and introducing the 
splitting function T> as in ((1.14(1 . we arrive in the relation 

Mx) + 0(M)]^ + + 004 + O(m)^ + (z, L[D]> = 0, (4.1) 

where the quantities 0{(£) as well as the coefficients of the first order linear differential operator 
L depend on § u ' s . To prove the claim now, it suffices to solve the vector field conjugacy problem, 
which ensures the structural stability of the vector field xq = (j^, (u, ^), 0^ on £ a ,p,T,r- The same 

conjugacy problem, only without the variable z, was dealt with by Sauzin, JJj (who referred to this 
problem as finding the characteristic vector field) , see also (12j , ^S] . The only difference here is the 
presence of the quantity z. However, as there is no differentiation in z in the "unperturbed" vector 
field xq, the quantity z enters the conjugacy problem as a parameter, and hence the resolution of 
the conjugacy is solely based on the invertibility of the operator + (cj, -g-) on Qa,p,T- Thus, the 
proof that the equation ((4.1(1 can be conjugated to 1(1.19(1 reproduces the proof of Lemma 4.4 in ^5] 
verbatim; we skip it, referring the reader to the latter or in fact any of the three above-mentioned 
papers. 

As the rest of the claims of Theorem ^ have been shown earlier in section 1 to follow from the 
claim (ii), this completes the proof of theorem. □ 



5. Proof of Theorem [21 

The proof follows the lines of the proof of the KAM theorem for semi-infinite cylinders in (15j . 
incorporating the dependencies in the "hyperbolic" variable z and can be in a sense compared to 
Graff's proof of the KAM theorem for whiskered tori, see e.g. J0]> |20j . 
Consider the differentiation operators 

A, = (u, Dp), D X(hU1 = X D S + D u . (5.1) 

The standard KAM theory depends on solvability of linear PDEs with the operator D u , in ^] the 
operator D\ LU was dealt with. 

Consider a perturbation of as in ((3.11(1 . The principal step in proving the structural stability 
of Hamiltonian the ((3.9|) is establishing the fact that the Hamiltonian if w is stable infinitesimally. 
This is done by solving the homological equation in the functional linearization of the problem 
(i.e. vindicating an "iterative lemma"). The standard Newton's iteration follows, see ^!5]> HO]- 
Parameter dependencies and smallness conditions were worked out for the case m = in the 
case m > makes no difference in this respect. Indeed, the estimates in the series of propositions in 
the Appendix, dealing with inversion of the first order differential operators involved, are in essence 
the same as they are in the latter paper, where the reader is referred for technical detail. 

The unknown quantities (S, b) appearing in ((3.10(1 exist, provided that one can solve the following 
set of equations (with the norm of the solution not exceeding the norm of the right-hand side by a 
factor, polynomial in the key parameters, such as analyticity loss): 
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[D Xo!U} + (z,AD z )]S = -f-(io,0+c, 

(5.2) 

[D Xo ,u + (z,AD z )]b = g + D*^ u (p,q)\j^(dS + i)+B(<p,x,z)b-\o-A!)Z. 

The system l)5.2|) arises by direct substitution of (|3.1U|) into ()3.9j) and omitting terms which are 

O2 ( I + |b| + |V"|). As far as the notation is concerned, £ = (£,0,0), Ao = (0, Ao,0) are n + l + m 

A y 

constant column-vectors and A = (0,0, Aq ) is a constant (n + 1 + tn) x m matrix. 

The role of the parameters c, Ao,Ao (in addition to £) is to ensure solvability of (|5.2|) within 
the framework of propositions in the Appendix, i.e. to guarantee that the right hand side is in 
the complement to the kernel of the operator D\ U} + (z,AD z ) on 25p(C) for the first equation 

and Qji"' A ((E) for the second one. Equivalently, after the canonical transformation \l/(a, £), the 
momentum-linear part of the Hamiltonian H w would acquire a term 

H UJ = c + \ e + {z,A z). (5.3) 

The term B(x, z)b can be described as follows. lib = (/3, b, b), describing the transformation of the 
(ip, x, z) variables respectively, then B(x,z)b contributes to the equation for the quantity b only, 
where it results in the term 

A T (s)b + bD s A T (s)z (5.4) 

in the right hand side. 

In order to solve the first equation in (|5.2|) . c is to be found, depending on the still unknown £, 
such that the right hand side, call it v^ o , have zero </j-mean (v§ ) = 0, cf. ()3.5j) . Recall that the 
mean it is taken by setting (s, z) = (—00, 0). 

No matter what £, such c clearly exists, so we can assume that the right hand side of the first 
equation has zero mean. Then Sq exists, in any space 23 p /(£), with p' < p, by Proposition IA.3I 
Observe that Sq is independent of £. 

The second equation in (|5.2j) comprises three (systems of) equations: for the quantities f3, b and 
b. First one considers the equation for (3 and finds £ such that the right-hand side, call it v&, has 
zero 99-mean, i.e. (vz) = 0. Note that the last three terms in the second equation in 1)5. 2|) do not 
appear in the equation for f3. Hence by the non-degeneracy assumption, 

i = -(DlH^{ P ,q) p=Q )- l Vp, 

where embraces the first n components of the n + l + m vector g + D^ rp H u} (p, q)\ p =odSo, member 
of the space S Q3p?' A '"^(£), for any p' < p. This also determines the constant c in ([5.3)1 . 

Furthermore, the (scalar) 6-component of the second equation in 1)5.2)1 is resolved by Proposition 
IA.5I The equation is not soluble without the condition 1)1. 9)1 . (The term constant Aoe is the only 
thing here to be added to Hamiltonian H w , because under condition ()1.9)) constants exhaust the 
kernel of the operator D\ 0:LU + (z,AD z ) on the space 2$p (<£)). 

Eventually, the equation for the quantity b is solved. This equation deserves special attention, 
so let us write it down explicitly as follows: 

[D\ 0iUJ + (z,AD z ) — A T ]b = Vy (5.5) 
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Let 

v~ b (cp,s,z) = v~ bfi ((p,z) +v^ 1 (ip,s,z), 

b((p,s,z) = b (ip, z) + bi(</>, s,z), 
where the quantity vr Jtp, — oo, z) = 0, and so it satisfies the estimate (|3.8|) . Therefore, the quantity 
bi(<p, s, z) G [2$*, (£)] m , corresponding to the right-hand side 1 exists, by Proposition IA.4I 
It remains to determine \)o(ip,z). Let 

\>o((p, z) = t> ,o(</>) + {z, bo,i(v?)) + 02(2; y), 

do the same expansion for the right-hand side Q . Then bo,o is found by Proposition IA. 21 cf. (|A.6|) . 

As for the term bo, 1 (</?), it is easy to see that the quantity {z, const.) is in the kernel of the operator 
in square brackets in 1)5.5)) . Hence the quantity Ao is introduced to ensure that the right hand side 

do not contain a constant multiple of z. This having been done, for all z, the right hand side 
(z, Q where Q ^(p) has zero mean, can be resolved by Proposition IA.11 (i). 

Finally, the component v^ Q2 = 02(z;ip) in the right-hand side of equation 1)5. 5 jl gets taken 

care of as follows. Consider a monomial z^ 1 . . . z^Uk^tp), with \k\ = k\ + . . . + k m > 2. Under the 
action of the operator (z,AqD z ) — Aq , taking into account the fact that Ao is diagonal, one gets 
some z^ 1 . . . z^h.QUk(ip), where the constant matrix Ao is diagonal and by the condition 1)1 .9|) is 
such that the real part of each diagonal entry is strictly positive, bounded away from zero uniformly 
in (k\, . . . , k m ) by some A > 0, which may be set equal to, say one tenth of the infimum in the 
right-hand side of 1)1.9)1 . Then the equation gets resolved term by term in the same way as is ()A.3)) 
in Proposition IA.21 the bound for the norm being uniform for all powers of z. 

The proof of infinitesimal stability of Hamiltonian H u will be complete after diagonalizing the 
constant matrix Ao + Ao by the linear transformation z — > Lz, where L is a constant near-identity 
matrix, such that L _1 (Ao + Ao)L is diagonal. This is possible as long as Ao is small enough. 

As we have mentioned earlier, this suffices to prove Theorem |^1 as one can now switch on the 
Newton's iteration procedure, see JH]) |20| . For the estimates, which would result in the forthcoming 
qualitative version of the theorem, with the smallness condition and parameter dependencies, see 

na. □ 

Quantitative statement of Theorem [2] 

We now present a quantitative statement of Theorem[21 The qualitative assumptions naturally look 
somewhat tighter than as stated in Theorem |^1 

Assumption 2' Suppose 3p = (a,T,p,r) > 0, as well as {pL,v) : < fj, < v <C 1, such that 
H u G Q3 K>P (T*£) and in the perturbation one has 

f G »p(£), g G 58^ A ' m) (£), |/| p < /i, \g\ p < ^ . (5.6) 
Regarding the terms in the expression h'-i. !J\) for H^, suppose 
i. <jj G M n satisfies JO) : 

ii. 3 A > 0, such that Vs G Rt,p the eigenvalues Xi(s), . . . , A m (s) of A(s) = Ao + Ai(s) satisfy 
A < mm(RXj(s)) < max(KAj(s)) < Aq — A, j = 1, . . . m; 
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Hi. Aq = diag(Xo t i, . . . , Ao, m) ) and Vk £ TH A 



^kjXoj - A 

3=1 



> A; 



iv. for any constant m x m matrix A, with ||A|| < A, the matrix Aq + A is diagonalizable; 

v. 3R,M>0, suchthaty(p,q)eM^+ 1 + m x(t p , \\ (D 2 U 2 (p; q))' 1 \\ < RT 1 , {{D^ip; q)\\ < M. 

Let < p' < p. Further without loss of generality assume that the quantities 5 = a — a', A = |p — 
p^A^M" 1 , H" 1 < 1. Theorem^ now vindicates the existence of a canonical transformation ^ 
such that {H UJ + V)o^ = H'^, where H' w 6 Q3 K / ;) y(T*(£) satisfies Assumptionl^lwith slightly modified 
parameters A' , A'(s), R' , M' . The quantitative results and parameter relations, cf. Assumption 
can be summarized as follows. 

Theorem 2' Under Assumption\$\ let k! = k/2 and 

? = inf( 7 (5 r " , A), r? = flinf(Af _1 ?A,i/). (5.7) 
There exists a constant C, depending only on ip, as well as the quantities n,T n ,ifj,p, k, such that if 



li < C~ 2 r, 2 



< 



(R/M) 2 A 2 [inf v)f 



(5.8) 



the following estimates hold: 
l§lp' 



< 



l b lp' 



AoVo-Aol < CfiinX)- 1 , A j|A^.-A 0li | < C^(r,\)-\ 
R- l \R'-R\ < CuirjqA)- 1 , M^M' - M\ < C/x^a) -1 . 



(5.9) 



The smallness condition q5.8|) is essentially the same as it was in |15j . It is determined by the 
estimates in the series of propositions in the Appendix only. These estimates coincide with the 
estimates in the corresponding propositions in the latter reference, where the resulting estimate, 
the analog of 1)5. 8|) is discussed in detail. 



6. Appendix 

The appendix contains a series of propositions necessary to resolve the infinitesimal conjugacy 
problem in the proof of Theorem [2] in this paper. The corresponding first order linear PDEs are 
solved by the method of characteristics; the proofs bear a close relation to lemmata in Chapter 5 
in Zehnder's work |2L)| . where the reader is directed for extra detail. 

The first result is adopted from [Ej- It is based on the classical result regarding the operator 
Du, due to Riissmann, ^Sj. The frequency uj is assumed to satisfy 1)1. 2j) . although this assumption 
is necessary only in the context of the operator D^. 

Proposition A.l i. For a function v £ *B a (T n ) with (v) = 0, the solution of the equation 
D^u = v exists in the space s 3 a /(T n ) for any a' < a. If a — a 1 = 5, ? = ^8 Tn , then 

Ha' ^ ?~ 1 b| (T . 
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ii. Let p = (a,T,p) and v G <B p (C), with (v) = 0. The solution of the equation D\ o ul u = v exists 
in QSp'(C) for any p' = (a',T,p) with < a' < a. If a — a' = S, q = inf (jS Tn , Aq 1 ), then 

in. For v G 23 p (6) ; there exists a real constant c, \c\ < \v\p, such that the solution of the equation 
D\ ,oj u = v — c exists in 5Sp/(C) and for the same ? as in (ii) one has 

Hp' < -^Hp- 
Proposition A. 2 Let p = (o~,T, p) and v G [53 p (C)] m . Consider the equation 

[D Xo ,u ~ Hs)]u = v, (A.l) 

where the matrix A(s) G [25 p (C)] m is such that any eigenvalue of the constant diagonal matrix 
Aq = A(— oo) = diag(\\, . . . , X m ) satisfies < cA < KAj < CA < Xq — cA for some c, C > 0. 
The solution of \A . 1\) exists in 23p(C) and with ? = inf(7<!> Tn , A), 

Hp < r l \v\ 9 . (A.2) 



Proof. The characteristic flow of D\ 0)LU is <f>t{y, s) = (ip + u>t, s + Xot), which clearly maps C p into 
itself. 

Decompose v(tp,s) = v (tp) +v 1 (p,s) (with v G [<B CT (T n )] m and v 1 G [<B*(e)] m ) and A(a) = 
Ao + Ai(s) in the sense of (|3.4|) . Seek the solution u(ip,s) = uo((p) + ui(p,s), expecting to find 
n G [*B CT (T n )] m and u x G [53j(C)] m . Then for u we have 

(D u - A )u = vo, (A.3) 

while u\ should satisfy 

D t ui - Ai(<^(s))ui = uiOt(</>> ( A - 4 ) 

where -Dj means differentiation along characteristics. The solution of equation (|A.3|) involves no 
small divisors and exists as long as the matrix Ao is non-singular and diagonalizable. It is assumed 
that Ao is diagonal, so the system of equations (|A.3|) gets separated into m equations: 

(D u - Xj)uoj(if) = v ,j(tp), j = l,...,m. (A.5) 

This results in an obvious bound |uo| CT ^ A _1 |uo| CT ) as each individual equation in (|A.5|) gets solved 
as the Fourier series with the coefficients 

Uoj, k = r J 3 !" r , fcGZ". (A.6) 
— Aj + t{k, u) 

Note that the constants c, C get absorbed into < symbols. 

For equation ()A,4|) let h(ip, s,t,t') solve the homogeneous equation 

D t h(ip,s,t,t') = AiOiO, s))h(<p, s,t,t'), 
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with h(ip,s,t,t) = 1. As Ai(4> t ((p, s)) = A\(s + Xot), one concludes that h does not depend on (p 
and moreover h(s, t, t') = h(s + Xot, s + Aq^). 

Moreover, for t' < t < one has the growth condition 

\h(8,t,H)\<e ox <t-*\ (A.7) 

Then, as (ip + ut',s + X t')\ < e s+Ao *'|ui| p! by definition of the space Q5p(C), cf. ()3.6j) . the integral 
in the right hand side of the representation 



u\(ip,s,t) = / h(s, t, t')vi(ip + cot', s + Xot')d£ (Ai 

J — oo 



converges absolutely for all t > 0, uniformly in s, with the bound |ui| p < (Ao — CA) 1 || - ui|| p , and «i 
will be a member of the space [53p(C)] m as is v\. □ 

Proposition A. 3 Let p = {cr,T,p,r) and v G 53 p (£), wrai/i (i>) = 0, let A(s) 6e sitc/i i/ia£ for all s, 
all its eigenvalues have positive real parts, bounded from zero by X > 0. The solution of the equation 

[D Xo> u + (z,AD z )]u = v (A.9) 

exists in 53 p /(£) for any p' = (a' ,T, p,r), with the bound hA.ty of Proposition VK~Jk 

Proof. The characteristic flow of the operator in square brackets in (jA.9|) is (f)t((p,s,z) = (ip + 
tot, s + Xot, ((z, s, t)), where C,(z, s, 0) = z and ( = A T (s + Xot)(. 

By positivity of Ao and the assumption on the eigenvalues of A, bounded in terms of A, the 
characteristic flow <j>t is well defined on (— oo,0] x (£ p , and we have an estimate 

IC(Ol<e" AM,| lC(t)l, i'<i<0. (A.10) 

After the decomposition v = Vo(cp, s) + (z, v\{ip, s, z)) and the same for u, the quantity uq is found 
after Proposition IA.11 

Furthermore, v\ G [<B p (£)] m still satisfies |t>i| p < \v\f p (recall that by convention r is absorbed 
into < symbols), besides 

{D XotlJJ + {z,KD z ))u 1 + Ku l = v l . (A.ll) 
Now let h(ip, s,t,t') solve the homogeneous equation 

D t h((p, s, z, t, t') = -A(4> t ((p, s))h(ip, s, z, t, t'), 

with h((p, s, z,t,t) = 1, where Dt is differentiation along characteristics. Clearly for t! < t < one 
has 

|%>, S ,2,M')l<e A ^\ (A.12) 

cf. (fA~7l) . 

Hence one can let 

ui(f,s,z)= h(cp,s,z,0,t)vi((f)t(cp,s,z))dt, (A.13) 



which guarantees that u\{ip, s, z) p < A 1 ||vi|| p as well as the fact that u\ G [53 p (£)] m . □ 
The following Proposition follows immediately from Propositions IA.21 and lA~3l 
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Proposition A. 4 Let p = (o~,T,p,r) and v £ [53p(£)] m . Consider the equation 

[Dx ,u + (z, A(s)D z ) - A T (s)]u = v, (A.14) 

where the matrix A(s) G [5S p (C)] m is such that for all s all its eigenvalues have positive real parts, 
bounded from above by CX < Xq — X, for some A > 0. 
The solution of £A.14 ) exists in [QSj((£)] m , with 

M P < A^Mp. (A.15) 

Proof: The characteristic flow (pt(f,s,z) = {ip + u>t,s + Xot, ((z, s, t)), with ((z, s,0) = z, of the 
operator Dt clearly maps £ p into itself. The solution of the homogeneous equation h(tp, s, z,t,t') 
satisfies estimate (|A.7|) . so by (|3.8|) it becomes possible to define 

u(ip,s,z,t)= / h(if, s, z, t, t')v [if + ojt' , s + Xot', C(z, s, t'))dt', (A. 16) 

J — oo 

which satisfies ()A.2|) . □ 

Proposition A. 5 Let p = (a,T, p,r), consider equation &A.91) . with v £ 53p(£). Suppose A(s) 6 
[53p(C)] m is such that all its eigenvalues have positive real part, bounded away from zero by some 
A > 0. In addition, suppose A(— oo) = Ao = diag{X\, . . . , A m ) and for any k £ Zip one has 

in 

\X - ^2 kjXj\ > X. (A.17) 

3=1 

There exists a real constant c, |c| < \v\p, such that the solution of \A.9\) exists in <Bp, (£) for any 
p' = (a',T,p,r) with the bound \A.fy) . 

Proof. The variation from Proposition IA.3( l (the characteristic flow is the same, Dt standing for 
differentiation along characteristics) is clearly the fact that 

v(tp,s,z) = V0 ^\ Z) +vx(fp,8,z), (A.18) 
X{s) 

So u also has to have a term u °^ s '^ ■ Substituting this term into (|A,9|) we get 

(-X D\nx(s + t) + D t )u {(t>t(<p,z)) = — — — -v (4> t (f, z)). 



Note that one can represent -Dlnx(s) = 1 + x( s ) w ( s ), with w(s) £ <B P (C), so the problem will 
reduce to Proposition IA.3I if we can solve the equation 

(-Xo + D u + {z,AD z ))uoWp,*)) = MM t P>*))- (A.19) 

Try uq as a monomial uq^,...,^^)^ 1 ••• z ^Ti with k £ Z™, substitute it in the latter equation, 
with the monomial fo,fe 1 ,...,fc m ( ( / 3 ) z i 1 • • • z m l 1Xi the right-hand side. This yields 
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Du + ^2 k j^j - ^0 UQ >kl> ..., km (<p) = Uo,fc 1 ,...,fc m (v), 
3=1 J 

which implies \uo\p < A _1 |t>o| p , by ()A.17|) . cf. (|A.6|) . 

The equation for u% with the right-hand side v\ from (|A.18|) is now amenable to Proposition IA.31 
the right hand side being v\ = v\ + Xow(s)uo((p, z). In general (vi) ^ and should be compensated 
by the constant c; it is not difficult to show that in fact c = D 2 ip(0)(vo) + (vi) (see |15| . Proposition 
B.4). □ 
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